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The  two-step  (7j  and  Tf)  deformation  behaviour  of  Ni3Al-based  single  crystals 
was  modelled  under  the  framework  of  a  new  constitutive  model  proposed  by 
Y.S.  Choi,  D.M.  Dimiduk,  M.D.  Uchic,  et  al.  [Phil.  Mag.  87  1939  (2007)].  A  new 
set  of  formulations  and  criteria,  which  identify  thermally  reversible  and 
irreversible  components  of  the  constitutive  variables  and  define  the  relative 
significance  of  those  components,  was  developed  and  implemented  within  the 
new  constitutive  framework.  The  simulation  results  well  captured  the  general 
qualitative  trends  of  the  flow  behaviour  upon  re-straining  at  T2  after  pre-straining 
at  7j  for  both  TX>T2  and  TX<T2.  Modelling  results  suggested  that  the 
dislocation  substructures  generated  at  7j  need  to  be  treated  as  partially  or  fully 
transferable  to  plastic  flow  at  T2 ,  at  least  through  the  early  stage  of  re-straining, 
to  capture  all  major  pre-strain  effects.  In  particular,  the  large  strengthening  effect 
at  T2  for  even  a  few  percent  of  pre-strain  at  7j  was  obtainable  only  by  controlling 
the  availability  of  mobile  dislocations  and  sources  at  T2. 


1.  Introduction 

The  partially  reversible  flow  stress  upon  re-straining  at  T2  after  pre-straining  at  T , 
(T2f  7j)  is  one  of  the  key  thermomechanical  features  of  Ni3Al-based  intermetallic 
single  crystals.  This  peculiar  behaviour  is  significant  since  it  holds  a  clue  to 
understanding  the  dislocation  behaviour,  its  link  to  evolving  substructures  during 
plastic  deformation  and  how  these  constituents  contribute  to  the  macroscopically- 
observed  flow  stress.  Since  Davies  and  Stoloff  [1]  first  utilized  the  two-step  (71  and 
Tf)  deformation  experiment  for  Ni3Al  to  identify  the  nature  of  plastic  flow, 
numerous  other  experimental  studies  have  sought  clarification  of  the  mechanisms 
governing  two-step  deformation  behaviour  of  Ll2  intermetallics  [2-13],  Ezz  and 
Hirsch  experimentally  verified  that  a  Cottrell-Stokes  law  [14]  holds  for  the 
deformation  structure  of  Ni3(Al,Hf)B  single  crystals,  which  shows  that  irreversible 
dislocation  storage  processes  are  partially  responsible  for  strain  hardening  [15]. 
However,  that  same  study  serves  to  illustrate  the  remarkable  and  exclusively 
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thermally  reversible  component  of  the  flow  stress.  This  implies  that  some  of  the 
stress-controlling  dislocation  configurations  are  thermally  reversible,  others  are 
irreversible  and  both  may  contribute  to  plastic  yielding,  sustaining  the  flow  stress 
and  setting  the  strain-hardening  rate  (SHR). 

Ezz  and  Hirsch  [12]  also  studied  the  effect  of  RT  pre-strain  ( 7)  <  T2)  on  the  flow- 
stress  anomaly  in  a  Ni3(Al,Hf,B)  single  crystal.  The  study  showed  that  a  pre-strain  at 
T\  significantly  changes  the  temperature  dependence  of  the  flow  stress,  as 
schematically  illustrated  in  figure  1  (identical  to  figure  6  in  [12]).  Here,  the 
temperature  dependence  of  the  yield  stress  for  the  RT-pre-strained  case  was  divided 
into  two  distinct  regimes:  ‘A’  (RT<7’<7’C)  and  ‘B’  (T>TC)  by  defining  a 
temperature  Tc,  at  which  the  yield  stress  ( rvp(T )),  upon  re-straining  at  T2  after 
pre-straining  at  RT,  is  equal  to  the  yield  stress  (r y(T))  without  RT  pre-straining.  In 
figure  1,  regime  A  exhibits  a  significantly  diminished  positive  temperature 
dependence  of  the  yield  stress  due  to  RT  pre-strain.  However,  in  regime  B  ryp(T) 
is  just  slightly  lower  than  ry(T).  For  compression  near  the  [0  0  1]  orientation,  Tc  was 
in  the  range  673-803  K,  depending  upon  the  magnitude  of  RT  pre-strain  [12]. 

In  contrast  to  the  numerous  experimental  studies  on  the  two-step  deformation 
behaviour  of  Ni3Al  alloys  [1-13,  15],  only  a  few  studies  attempt  theoretical  modelling 
of  this  behaviour.  Both  Ezz  and  Hirsch  [12]  and  Greenberg  and  Ivanov  [16]  proposed 
theoretical  frameworks  for  the  two-step  deformation  behaviour.  However,  both  of 
these  frameworks  were  solely  qualitative  and  quantitative  modelling  has  never  been 
attempted. 

In  [17],  the  present  authors  proposed  a  comprehensive  crystallographic 
constitutive  model  for  Ni3Al-based  intermetallic  single  crystals.  The  present 
study  used  that  constitutive  framework  to  model  the  two-step  deformation 
behaviour,  viz.  the  thermal  reversibility  of  the  flow  behaviour,  of  Ni3Al-based 
single  crystals.  Here,  both  the  physical  basis  for  the  constitutive  model  and  the 
understanding  gleaned  from  its  framework  serve  as  the  conceptual  framework  for 
modelling  the  two-step  deformation  behaviour.  Those  concepts  are  summarized  in 


Figure  1 .  Schematic  plot  illustrating  the  effect  of  RT  pre-strain  on  the  positive  temperature- 
dependence  of  yield  stress  as  adopted  from  figure  6  in  [12].  The  plot  was  based  on  experimental 
observations  from  compression  of  near-[001]-oriented  Ni3(Al,1.5%  Hf)0.2%  B  single  crystals 
by  Ezz  and  Hirsch  [12], 
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section  2  of  this  report.  A  thorough  review  of  previously  published  experimental 
and  theoretical  work  on  the  two-step  deformation  behaviour  of  Ll2  single  crystals 
formed  the  basis  for  this  study  (refer  to  [1-13,  15,  16]).  Based  upon  this  review, 
key  concepts  for  the  thermal  reversibility  of  Ni3Al-based  single  crystals  were 
identified  and  used  to  represent  the  effect  of  pre-strain  within  the  framework  of 
our  constitutive  model.  An  interpretation  of  those  concepts  for  the  purposes  of 
implementing  the  present  model  is  discussed  in  section  3.  Finally,  numerical 
simulation  was  performed  to  investigate  an  incorporation  of  the  pre-strain  effect 
in  the  constitutive  model.  The  results  are  discussed  within  the  context  of  their 
physical  implications  in  sections  4  and  5. 


2.  Constitutive  model  proposed  in  part  I 

A  summary  of  the  new  constitutive  descriptions  for  Ni3Al-based  single  crystals  can 
be  found  in  section  3.3  of  [17].  That  constitutive  model  was  structured  upon  an 
understanding  of  the  extensive  literature  describing  physically  verified  micromechan¬ 
isms  associated  with  the  anomalous  flow  behaviour  of  Ni3Al-based  Ll2  inter- 
metallics.  For  those  materials,  upon  plastic  straining,  mobile  screw-character 
superdislocations  are  ‘effectively’  instantaneously  exhausted  by  forming  Kear- 
Wilsdorf  locks  (KWLs)  [18,  19].  The  locks  result  from  intrinsic  elastic  and  fault- 
energy  related  forces,  driving  cross-slip  of  screw-character  superdislocations  from  the 
octahedral  plane  to  the  cubic  plane,  resulting  in  a  range  of  superlattice-partial 
dislocation  separation  distances  in  the  cubic  plane.  The  locking  mechanism  also 
results  in  a  spectrum  of  non-screw-character  superdislocation  segments,  called 
macrokinks,  remaining  glissile  on  the  octahedral  plane.  Thus,  elementary  glide  loops 
adopt  an  anisotropic,  three-dimensional  nature  while  maintaining  effectively 
octahedral  glide.  Such  exhaustion  limits  the  availability  of  mobile  dislocations, 
which  results  in  an  STIR  that  well  exceeds  any  rate  expected  from  ordinary  FCC 
single  crystals.  The  exhaustion  of  mobile  dislocations  and  its  influence  on  the  SHR 
and  flow-stress  anomaly  was  theoretically  treated  by  Caillard  [20]  and  numerically 
modelled  by  the  present  authors  [21].  The  recent  slip-trace  analysis  of  Ni3Al  using 
atomic  force  microscopy  by  Fikar  et  al.  [22]  and  Coupeau  and  Bonneville  [23]  also 
confirmed  the  importance  of  the  mobile-dislocation  exhaustion  to  the  SHR  and 
flow-stress  anomalies. 

The  dislocation  substructure  associated  with  this  exhaustion  mechanism  is 
typified  by  long,  straight  locked  screw-character  superdislocations  (KWLs) 
connected  by  a  series  of  macro-kinks  (MKs)  and  cross-slip  jogs  (CJs)  on  the 
cubic -plane  steps  connecting  the  octahedral  planes.  The  new  constitutive  model 
represented  such  dislocation  substructures  using  two  important  constitutive  structure 
parameters,  the  critical  MK  height  (/?,.)  and  the  average  MK  height  (/?),  both 
determined  in  the  octahedral  plane._  The  parameter  hc  is  the  critical  height 
above  which  MKs  are  mobilized  and  h  is  the  average  distance  that  screw-character 
dislocations  glide  in  the  octahedral  plane  before  they  are  locked.  Once  mobile  screw- 
character  dislocations  are  ‘effectively’  immobilised  by  forming  KWLs,  further 
plastic  straining  is  facilitated  by  two  mechanistic  events:  the  mobilization  of  ‘weak’ 
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MKs  (edge  character)  and  self  (athermal)  unlocking  of  ‘weak’  KWLs  (screw 
character).  The  new  constitutive  model  treated  the  former  using  a_fraction  of 
mobilized  MKs  (/mK),  which  is  a  function  of  hc  and  h  (fmK  —  exp(— hc/h))  based  on 
the  statistical  treatment  of  MK-height  populations  and  the  latter  using  a  threshold 
stress  (r,),  which  is  the  stress  required  to  unlock  incomplete  KWLs.  When  the  applied 
stress  is  not  high  enough  to  reach  r h  viz.  in  the  regime  dominated  by  exhaustion  of 
screw-character  dislocations,  the  flow  stress  and  area  swept  by  dislocations  is  limited 
primarily  by  the  availability  of  mobile  MKs  (/mK)  (hence  the  mobile  MK  density 
pmK).  Here,  the  average  MK  height  (/;)  decreases  exponentially  with  increasing 
temperature  as  its  magnitude  is  inversely  proportional  to  the  frequency  of  cross-slip 
locking,  which  is  a  strong  positive  function  of  temperature.  This  means  that,  in  the 
anomalous  temperature  domain,  one  can  expect  a  higher  frequency  of  cross-slip 
locking  with  increasing  temperature,  leading  to  a  decrease  in  the  average  value  of 
MK-height,  which  eventually  limits  the  mobile  MK  density  (pniK)-  Clearly,  it  is 
important  to  accurately  represent  the  MK  behaviour  in  any  model  that  describes  the 
flow  stress  after  pre-straining  at  different  temperatures. 

The  new  constitutive  model  for  Ni3Al-based  single  crystals  is  distinct  from 
conventional  constitutive  models  for  fee  crystals.  A  major  difference  is  the 
minimal  contribution  of  forest-dislocation  obstacles  in  the  new  model  (at  least 
within  the  small  strain  regime).  In  other  words,  the  area  swept  by  mobile 
dislocations  is  primarily  limited  by  ‘effectively’  trapping  screw-character  disloca¬ 
tions  via  KWLs  and,  thereby,  limiting  the  MK  height  distribution  and  mobile 
density.  For  this  regime,  the  mobile  density  is  so  limited  that  the  stress  must  rise 
sharply  to  maintain  an  ample  supply  of  dislocations  to  satisfy  the  strain  rate. 
Such  a  mechanism  is  in  stark  contrast  to  what  is  envisaged  for  fee  metals, 
wherein  the  area  swept  is  limited  by  dislocations  interacting  with  forest- 
dislocation  obstacles.  For  the  forest-hardening  case,  the  material  exhibits  a 
much  lower  stress  sensitivity  of  the  mobile  density;  thus,  there  is  an  abundance  of 
mobile  dislocations  at  comparatively  modest  stresses. 

This  KWL-formation,  MK-motion  mechanism  is  expected  to  dominate  in  the 
small-strain  (up  to  3-5%  of  the  axial  strain)  exhaustion-hardening  regime.  Thus,  in 
this  regime,  the  flow  stress  of  Ni3Al  may  not  be  quantified  through  the  observed 
dislocation  density  and  the  Taylor-type  relationship  known  to  hold  for  fee  metals. 
This  was  indeed  the  major  conclusion  drawn  by  Kruml  et  al.  [24]  from  their 
measurements  of  flow  stress  and  the  dislocation  density  for  Ni3Al  alloys.  The  new 
constitutive  model  treated  these  concepts  numerically  through  equations  describing 
the  dislocation  substructure  resulting  from  such  mechanistic  processes.  Again,  the 
two  constitutive  parameters,  hc  and  h,  play  key  roles  in  representing  the  ‘strength’ 
and  the  temperature  dependence  of  such  a  dislocation  substructure.  In  particular, 
these  two  parameters  are  critical  for  representing  the  effect  of  pre-straining  at  T |  (i.e. 
the  influence  of  the  7)  dislocation  substructure  on  the  T2  deformation)  when 
modelling  two-step  (7)  and  T2 )  deformation  behaviour,  the  details  of  which  are 
discussed  in  the  next  section.  Nonetheless,  having  drawn  focus  to  the  primary 
importance  and  parameters  describing  MK  behaviour,  one  cannot  ignore  fcc-like 
contributions  to  strain  hardening.  Forest-like  obstacles  pose  a  resistance  to  MK  glide 
and  contribute  to  irreversible  components  of  flow  stress  during  two-step  deformation 
sequences,  as  detailed  in  what  follows. 
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3.  Incorporating  the  pre-strain  effect  in  the  constitutive  model 

A  literature  survey  [1-13,  15,  16]  and  interpretation  of  those  reports  identified  key 
phenomenological  aspects  of  the  two-step  deformation  behaviour  of  Ni3Al-based 
single  crystals;  these  key  points  are  listed.  For  the  sake  of  clarity,  the  word  ‘original’ 
is  used  hereafter  to  refer  to  the  material  state  that  has  no  pre-strain  history. 

(i)  For  7b  >  7b  experiments,  the  magnitude  of  pre-strain  at  T  {  was  usually  very 
small,  less  than  about  2.5%  shear  strain.  Even  though  the  flow  stress 
measured  when  re-straining  at  T2  was  lower  than  that  for  original  straining 
at  T\,  the  small  magnitude  of  pre-strain  at  7)  gave  rise  to  a  significant  flow- 
stress  increment  at  T2.  This  is  expressed  by  Ar  =  r^-r^j,  where  and  t(0) 
are  the  T2  flow  stresses  with  and  without  pre-strain  at  7),  respectively.  The 
observed  flow-stress  increment  at  7b,  Ar,  ranged  from  45  to  170%  of  r(o)  [5, 
6,  9,  13]  depending  upon  the  loading  orientation,  values  of  7)  and  T2,  the 
magnitude  of  pre-strain  and  the  alloy  composition.  The  magnitude  of  the 
flow  stress  determined  from  strain  hardening  at  Tx  was  observed  to  be 
comparable  to  Ar  at  7b.  However,  from  only  a  few  percent  pre-strain  at  7) 
ordinary  fcc-like  hardening  alone  is  not  expected  to  generate  such  a  large 
strengthening  effect  (Ar)  at  7b.  This  implies  that  strain  hardening  at  7)  may 
involve  complex  dislocation  substructures  whose  strength  contributions 
originate  from  various  sources  beyond  ordinary  fcc-like  hardening. 
Furthermore,  these  substructure  contributions  may  be  partially  or  fully 
transferred  to  the  deformation  at  7b,  hence  raising  the  magnitude  of  At. 

(ii)  For  T\<T2  experiments,  two  distinctive  regimes  were  identified.  In  the  case 
that  a  difference  between  7)  and  T2  (AT—  T2—  7b)  is  relatively  small  (i.e. 
regime  A),  a  pronounced  flow-stress  increase  was  found  at  7b  from  RT  pre¬ 
strain  experiments  [12].  This  flow-stress  increment  at  7b  dramatically 
decreased  with  increasing  AT  (hence  7b),  resulting  in  a  diminished  positive- 
temperature  dependence  of  the  flow  stress  at  T2  [12].  In  the  case  of  a  large 
AT  (i.e.  regime  B),  the  flow  stress  at  7b  decreased  compared  to  the  original 
flow  stress  at  7b  by  10^10%,  thus  initiating  a  micro-yield  transition  typically 
of  a  magnitude  less  than  a  few  percent  in  the  small  strain  regime  [5,  13]. 
These  behaviours  are  sensitive  to  the  magnitude  of  pre-strain,  the  loading 
orientation,  AT  and  the  alloy  composition. 

(iii)  For  both  7'1>7b  and  TX<T2,  the  strain-hardening  rate  of  material  re¬ 
strained  at  T2  tended  to  follow  the  original  flow  behaviour  at  7b.  This  trend 
was  more  evident  for  the  case  of  T\  >  7b,  wherein  crystals  exhibited  almost 
identical  strain-hardening  rates  at  7b,  with  and  without  pre-strain,  after  an 
early  substructure-rearrangement  stage  [6,  9,  13].  However,  such  a  trend  was 
also  clear  for  the  case  of  7)  <  T2,  wherein  strain-hardening  rates  for  crystals 
re-strained  at  T2  tended  to  quickly  recover  the  original  strain-hardening 
rate  at  7b  after  a  short  re-adjustment  (i.e.  the  micro-yield  transition)  at  the 
early  stage  of  re-straining  at  7b  [5,  13].  In  this  sense,  the  overall  plastic  flow 
of  Ni3Al-based  single  crystals  can  be  described  as  a  thermally  reversible 
process,  the  degree  of  which  depends  upon  the  balance  between  the 
strengths  of  reversible  and  irreversible  contributions  to  plastic  flow. 
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The  new  constitutive  model  for  Ni3Al-based  single  crystals  was  structured  from 
equation  (3)  in  [17],  which  can  be  re-written  as: 

Y=Yk  +  Yw  =  bpmKvK  +  bpmwvw.  (1) 

One  can  refer  to  [17]  for  detailed  descriptions  of  each  parameter.  In  equation  (1), 
Yw  reflects  the  plastic  flow  sustained  by  screw-character  dislocations  typified  by 
cross-slip  locking  and  athermal  unlocking  of  KWLs  having  a  variable  cross-slip 
distance  in  the  cubic  plane,  w.  This  component  of  the  strain  rate  (stress)  can  be 
understood  as  thermally  reversible  flow  since  its  basic  nature  is  set  by  intrinsic 
driving  forces  and  the  ‘ effectively ’  instantaneous  rate  of  KWL  formation.  The 
locking  rate  sets  yw  and  it  reflects  the  present  temperature  in  a  fully  reversible 
fashion.  This  implicitly  means  that,  in  equation  (1),  yjc  is  the  term  influenced  by 
pre-strain  history  at  7)  under  two-step  (7)  and  T2)  deformation.  Thus,  the  plastic 
strain  rate  upon  re-straining  at  T2  after  pre-straining  at  7)  can  be  expressed  by: 

Yip)  —  YK(p)  +  Yw(p)  —  bpmK(p)VK(p )  +  Yw,  (2) 

where  all  parameters  having  subscript  (p )  are  in  a  state  at  T2  after  being  influenced  by 
pre-straining  at  7) .  Parameters  having  no  subscript  (p)  are  in  their  characteristic  state 
for  T2  regardless  of  pre-strain  history  at  7) .  Thus,  in  equation  (2),  the  influence  of 
pre-strain  on  plastic  flow  at  T2  is  reflected  through  pm K(/))  and  vK(p). 

Based  upon  the  viewpoint  represented  by  equation  (2)  and  point  (i)  at  the  start  of 
this  section,  we  envisaged  that  two  types  of  substructures  (forest  obstacles  and  MK 
distributions),  which  were  generated  and  responsible  for  strain  hardening  at  7), 
affect  plastic  flow  at  T2.  The  first,  ‘forest’  obstacles  developed  and  stored  at  7),  poses 
a  resistance  against  the  MK  glide.  That  glide  resistance  was  treated  through  z y 
(equation  (17)  in  [17]),  based  upon  a  re-interpretation  and  simplification  of  the 
concept  of  flow-stress  partitioning  proposed  by  Ezz  and  Hirsch  [12,  15,  25].  This 
irreversible  fcc-like  obstacle  storage  was  expressed  as  7/(7),  /i),  where  y\  is  the  shear 
pre-strain  at  7)  and  it  is  assumed  to  play  the  same  role  as  obstacles  against  the  MK 
glide  at  T2.  Following  this  assumption,  Vkq>)  in  equation  (2)  can  be  expressed  by: 


Bp)  rK(p)  u 

vm  =  B  b. 

(3) 

And 

TK(p)  —  T0  ”£/(/?)  5 

(4) 

rAp)  =  Ty(7’i,y1)  +  T/, 

(5) 

where  B  is  the  drag  coefficient,  and  r ^  and  re  are  the  RSS  at  T2  and  the  CRSS  on  the 
octahedral  plane,  respectively.  In  equation  (4),  xK  is  the  slip  resistance  against  the 
MK  glide.  The  evolution  of  z y  at  T2  in  equation  (5)  was  set  to  follow  a  fcc-like 
parabolic-hardening  law  as  described  in  equation  (18)  in  [17]. 

The  second  substructure  type  affecting  the  flow  behaviour  at  T2  is  the  MK 
height  configurations,  which  were  generated  at  7).  In  the  new  constitutive  model 
[17],  plastic  flow  was  controlled  by  a  limited  availability  of  mobile  dislocations 
from  both  screw-character  dislocations  being  ‘ effectively ’  instantaneously 
exhausted  into  KWLs  and  edge-character  dislocations  having  a  mobility 
controlled  by  an  available  MK  fraction  (see  section  3.3  in  [17]  for  details). 
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Here,  the  former,  namely  the  formation  of  KWLs,  was  treated  exclusively  as  a 
thermally  reversible  event  through  equations  (1)  and  (2)  and  the  associated 
equations  given  in  [17].  However,  the  dislocation  substructure  associated  with  the 
MK-height  configurations  may  partially  or  fully  remain  after  pre-straining  at  T j 
and  probably  influences  plastic  flow  at  T2.  In  other  words,  the  development  of 
new  MK-height  configurations  at  T2  may  be  influenced  by  the  built-in 
substructures  associated  with  MK-height  configurations  at  7)  and  yx.  Also, 
viewpoint  (iii)  stated  above  suggests  that  such  influence  of  the  built-in 
substructure  may  be  temporary,  ranging  from  the  onset  of  plastic  flow  until  the 
newly-generated  substructure  at  T2  gains  full  control  of  plastic  flow  over  the  old 
substructure.  All  of  these  effects  were  incorporated  into  the  current  model 
through  pmK(rt  in  equation  (2).  From  equations  (5)  and  (6)  in  [17],  praK(/))  can  be 
expressed  as 


PmK(p)  —  PKtot  exp 


(6) 


where  pKtot  is  the  total  MK  density.  As  mentioned  in  the  previous  section,  in  the  new 
constitutive  model  [17],  the  dislocation  substructure  associated  with  MK-height 
configurations  was  represented  by  two  key  parameters,  hc  and  h.  In  equation  (6), 
these  two  parameters  were  used  as  carriers  of  substructure  between  7)  and  T2  to 
account  for  the  effect  of  built-in  MK-height  configurations  at  7).  In  the  new 
constitutive  model,  both  hc  and  h  were  assumed  to  asymptotically  decay  with  plastic 
strain,  as  expressed  by  equations  (7)  and  (8)  in  [17].  Similarly,  the  evolution  of  hc(p) 
and  in  equation  (6)  were  expressed  by: 


dhc(p) 

®hc(p)  | 

(  hc(p) 

and 

d  Y 

Vlco(p)  - 

d%) 

~9h(p)  [ 

{  h(p)  —  ‘ 

d  y 

\ho(p)  — 

hs 


(7) 

(8) 


where  dhc(/>)  and  0fi{pj  are  initial  decay  rates  for  hc(p)  and  h^),  respectively,  and  hcs  and 
hs  are  the  final  saturation  values.  From  equations  (6)  through  (8),  one  needs  to 
implicitly  understand  that  the  influence  of  built-in  MK-height  configurations  was 
incorporated  in  modelling  through  the  adjustment  of  initial  values  of  critical  and 
average  MK  heights  (hco(p)  and  h0(P),  respectively),  and  their  evolution  (&hc(p)  and 
with  strain  at  T2.  The  exact  arithmetic  forms  (equations  (7)  and  (8)) 
representing  a  smooth  evolution  of  substructure  between  states  7)  and  T2  are  a 
constitutive  assumption.  The  assumption  was  chosen  because  it  offers  a  reasonably 
well-behaved  representation  of  the  understood  MK  phenomenology  within  the 
constitutive  model,  even  though  this  aspect  of  MK  behaviour  has  never  been 
quantified. 

Figure  2  schematically  illustrated  these  parameters  at  T2  for  the  cases  with  and 
without  pre-strain  at  7).  Figure  2  takes  the  evolution  of  h,  at  T2  as  an  example.  Note, 
however,  that  the  same  illustration  is  also  applicable  to  the  evolution  of  h  at  T2.  Due 
to  pre-strain  at  7),  the  initial  critical  MK  height  at  T2  was  assumed  to  decrease  to 
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Figure  2.  Schematic  plot  showing  the  evolution  of  hc  at  T2  for  cases  with  (two  solid  lines) 
and  without  (broken  line)  pre-strain  at  7).  Dark  solid  lines  indicate  h c(P)  evolution  with 
different  initial  decay  rates,  @hc  and  9 hc(p),  respectively. 


hCO(p )  from  hco.  In  the  same  way,  the  initial  average  MK  height  at  T2  is  expected  to 
decrease  to  h0<j,)  from  h„.  The  values  of  hco(p)  and  depend  upon  those  of  hc(T\,  y i) 
and  h(T\,y\),  which  are  the  critical  and  average  MK  heights,  respectively,  at  the  end 
of  pre-straining  at  T\.  Here,  to  determine  hCO(j,)  and  h0(p),  two  preliminary  parameters 
Hco  and  //„  were  defined  by: 

and  Hco  ~  hc{T\,y\) +flc[hc0{T\)  -  hc{T\,y\)\  (9) 

H0  =  h(T\ ,Yi)  +fha[ho  -  h(T\ ,Yi)l  (10) 

where  /hc  and  /ha  are  the  adjustable  parameters  having  a  range  from  1  to  0, 
determining  the  relative  degree  of  influence  of  the  MK-height  related  built-in 
substructure  at  7j  on  the  development  of  the  MK-height  substructure  at  T2.  They 
represent  the  potency  of  the  unknown  physical  effects  of  partially  relaxing  the  7j 
substructure  on  unloading  and  of  transitioning  to  a  new  substructure  upon  straining 
at  T2.  The  highest  degree  of  influence  is  expected  for  the  values  fhc  and  /ha  equal  to  0. 
Through  equations  (6)  (10),  the  present  modelling  adopted  the  concept  of  viewpoint 
(iii),  which  accounts  for  the  pre-strain  effect  being  a  temporary  effect  only  in  the 
early  stage  of  plastic  flow  at  T2.  This  can  be  incorporated  in  modelling  by  letting  hc(p) 
and  h^)  follow  the  same  strain  evolution  as  that  of  hc  and  h,  respectively,  after  the 
initial  re-adjustment  stage  of  plastic  flow  at  T2.  For  this,  the  initial  decay  rates  0hc(/,) 
and  6p(/j}  for  hc(p)  and  are  required  to  decrease  to  values  lower  than  0hc  and  0 
respectively,  when  hcofp)  and  h(,(p)  are  lower  than  hco  and  h0,  respectively.  This  was 
schematically  illustrated  in  figure  2  by  a  dark  solid  line  showing  the  evolution  of  hc(p) 
having  0hc(p)  lower  than  0hc. 

For  the  determination  of  hco^  and  h,,^  and  their  corresponding  initial  decay 
rates  at  T2  the  following  equations  and  criteria  were  used: 

hco<j>)  —  Hco  and  9hr(p)  —  (1  .foc'j&hc?  if  Hco<^hco  (11) 

hcoip)  —  hco  and  0/,f(p)  =  6hc,  otherwise 
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and 

^o(p)  —  o  and 

ho(P)  —  h0  and 

Here,  hco  and  h„  are  the  original  critical  and  average  MK  heights,  respectively,  at 
the  onset  of  plastic  flow  at  T2.  These  were  defined  by  equations  (9)  and  (11), 
respectively,  in  [17].  In  equations  (11)  and  (12),  f9c  and  f9a  are  another  set  of 
adjustable  parameters  ranging  from  0  to  1L  depending  on  the  differences  between  hco 
and  ltCO(p)  {J\hco  —  hco-  hco(pf)  and  between  ha  and  h„(P)  (A h0  —  h0  —  h0(P)),  respectively. 
In  general,  as  A hco  and  A h0  increase,  the  corresponding  f9c  and  f9a  are  expected  to 
increase  to  ensure  a  short  re-adjustment  period  in  the  early  stage  of  plastic  flow  at  T2. 
Note  that  equations  (11)  and  (12)  are  effective  only  when  Hco  and  H0  are  smaller 
than  hco  and  h0,  respectively,  since  otherwise  the  original  hco  and  ha  already  take  full 
control  of  plastic  flow  even  at  the  initial  stage  of  re-straining  at  T2.  From  equations 
(9)  through  (12),  one  should  notice  that  /hc  and  /ha  control  the  degree  of  thermal 
irreversibility,  while  f9c  and  f9a  control  the  duration  of  the  influence  of  the  built-in  T \ 
substructure  at  T2.  The  values  of  /hc  and  /ha  being  close  to  1  bias  flow  towards  the 
reversible  nature  of  the  substructure  associated  with  MK-height  configurations  and 
vice  versa. 


%)  =  0  ~foa)6h,  if  H0<K 
%)  =  0h>  otherwise. 


4.  Simulation  outline 

The  two-step  (7)  and  73)  deformation  behaviour  of  Ni3Al-based  single  crystals  was 
simulated  using  the  formulations  described  in  the  previous  section.  These 
formulations  were  incorporated  into  the  framework  of  the  new  constitutive  model 
proposed  in  [17],  and  7)  and  T2  compression  simulations  were  performed  for  [1  2  3]- 
oriented  Ni3(Al,0.25%  Hf)  single  crystals  in  an  anomalous  temperature  regime  from 
300  to  900  K.  For  the  case  of  7)  >  73,  the  [I  2  3]  compression  was  simulated  up  to  1% 
axial  pre-strain  at  773  K  and  the  subsequent  re-straining  simulation  was  performed 
at  300  K.  For  the  case  of  7)  <  73,  however,  the  [1  2  3]-compression  simulation  was 
performed  up  to  the  axial  pre-strain  levels  of  5  and  10%  at  300  K  and  the  subsequent 
re-straining  simulations  were  carried  out  at  temperatures  from  300  to  900  K.  All 
simulations  were  preformed  at  a  strain  rate  of  5  x  10“5/s.  In  the  present  study,  four 
new  variable  parameters  were  introduced  in  the  previous  section,  /hc,  /ha,  f9c  and  f9a, 
and  used  for  parametric  studies.  The  values  of  all  other  input  parameters  were  the 
same  as  those  used  in  [17]  (see  tables  1  and  2  in  [17]  for  details).  A  preliminary 
parametric  _study  showed  that  f9a  expressed  as  an  asymptotic  function  of  A ha 
<f9a  =  0.5[Aho/20(A)]2/3)  delivered  a  reasonable  qualitative  dependence  off9a  on  A h„ 
with  varying  73  for  the  simulations  of  7)  <  73.  Note  that  the  same  numerical 
framework  for  the  two-step  (7)  and  73)  deformation  (equations  (2)  ( 1 2))  was  used 
for  both  cases  of  T\>T2  and  T\<T2.  The  simulation  results  are  discussed  based 
upon  qualitative  comparisons  with  the  experimental  data  and  the  physical 
metallurgical  implications  of  those  results  were  interpreted  in  the  context  of  previous 
anomalous  flow  theories. 
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5.  Simulation  results  and  discussion 
The  case  of  T/>  T2 

Figure  3  shows  the  simulated  [1  2  3]-compressive  flow  curves,  which  represent 
re-straining  at  300  K  after  being  pre-strained  by  1%  axial  strain  at  773  K,  for 
different  values  of /ha  from  0.1  to  0.4.  The  original  flow  curves  simulated  at  300  and 
773  K  are  also  shown  in  figure  3  for  comparison.  Values  of  fhc,fec  and  fga  were  set  to 
0,  0.1  and  0.5,  respectively,  for  the  current  simulations.  Flow  curves  simulated  for 
re-straining  at  300  K  well  represented  the  general  trend  of  the  partially  reversible 
flow  for  this  class  of  alloys  [5,  9,  13].  The  magnitude  of  i/(7’1,y1)  (equation  (5)) 
obtained  from  a  1%  axial  pre-strain  at  773  K  was  just  a  few  MPa,  which 
was  insufficient  for  obtaining  a  large  flow  stress  increment  (Ar=rw-r(„))  at  300  K, 
as  was  experimentally  observed.  The  data  from  Shi,  et  al.  [13],  revealed  an  almost 
100%  increase  of  the  flow  stress  for  experimental  tests  similar  to  these  simulations. 
In  the  present  simulations,  such  a  large  value  of  Ar  turned  out  to  be  achievable  only 
by  controlling  the  influence  of  built-in  7j  substructures,  specifically,  through 
decreasing  the  magnitude  of /ha  in  equation  (10).  The  low  magnitude  of  /ha  serves  to 
enhance  the  reduction  of  the  average  MK  height  (]i(P))  at  T2  (due  to  pre-straining  at 
7j),  which  ends  in  a  decrease  of  available  mobile  MK  density  (pmK(/>))  at  73  (through 
equation  (6)).  By  changing  the  value  of  /ha  to  0.4,  0.3  and  0.2  in  turn,  Ar  was 
increased  to  approximately  48,  77  and  130%  of  the  original  flow  stress  at  300  K, 
respectively.  These  simulations  strongly  suggest  that  a  large  Ar  for  the 
case  of  7j  >  T2  arises  mainly  from  the  reduced  availability  of  mobile 


Re- 

strained 
at  300  K; 

fha 


Figure  3.  Simulated  axial  stress  vs.  axial  strain  curves  for  different  values  of /hli  from  0  to  0.4 
re-strained  at  300K  after  being  pre-strained  by  1%  axial  strain  at  773  K.  Values  of  fha  fee  and 
fea  were  set  to  0,  0.1  and  0.5,  respectively,  for  the  current  simulations.  Original  axial  stress  vs. 
axial  strain  curves  simulated  at  773  and  300  K,  respectively,  are  also  shown  for  comparison. 
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dislocations  at  T2.  Qualitatively  speaking,  the  mobile  dislocation  density  in  the  early 
stage  of  plastic  flow  is  reduced  by  the  influence  of  the  relatively  ‘strong’  MK-height 
configurations  that  were  built  during  pre-straining  at  7).  This  effect  was  treated 
by  pmK(p)  of  equation  (6)  through  hc(p)  and  h(p)  (equations  (7)  and  (8))  in  the  current 
framework  of  modelling.  The  simulation  results  also  partly  agree  with  the  idea 
proposed  by  Shi  et  al.  [13],  wherein  they  postulated  that  a  significant  reduction 
of  dislocation  sources  occurs  at  T2  as  the  result  of  higher  temperature 
pre-straining  at  7). 


The  case  of  Tj  <T2 

Figure  4  shows  the  simulated  0.2%-offset  shear  stress  (r0.2%)  vs.  temperature  (T) 
plots  for  the  cases  pre-strained  by  5  and  10%  at  300  K,  followed  by  re-straining  at 
temperatures  from  300  to  900  K  (this  will  be  referred  to  as  ‘pre-strained’  hereafter). 
Figures  4a  and  b  are  the  simulation  results  for  two  different  values  of /ha,  0.1  and  0, 
respectively,  where  applicable.  The  simulated  original  r0.2%-T  plot  (figure  2  in  [17]) 
was  also  included  for  comparison.  Flere,  values  of  /hc  and  f6c  were  chosen  to  be  0.2 
and  0.4,  respectively,  for  the  case  of  5%  pre-strain,  and  0.45  and  0.5,  respectively,  for 
the  case  of  10%  pre-strain,  and  f6a  was  set  to  asymptotically  vary  with  A h0  using  a 
formulation, /e„  =  0.5[A/jo/20(A)]2/3,  where  applicable.  The  simulation  result  well 
captured  the  general  trend  for  the  effect  of  RT  pre-strain  on  the  positive  temperature 
dependence  of  the  flow  stress  (schematically  illustrated  in  figure  1),  based  upon  the 
experimental  work  by  Ezz  and  Hirsch  [12].  A  significantly  reduced 
positive  temperature  dependence  of  rQ.2%  was  obtained  in  simulated  r0.2 %—T  plots 
(figures  4a  and  b)  in  the  temperature  regime  below  Tc  (regime  A),  where  Tc  was 
about  600  and  650  K  for  the  cases  of  5  and  10%  pre-strain,  respectively.  The  major 
difference  in  the  simulated  r0.2%-T  plots  of  figures  4a  and  b  is  the  almost  complete 
absence  of  the  positive  temperature  dependence  of  r0.2%  in  regime  A  (particularly  for 
10%  pre-strain)  for  the  case  of  zero /ha  (figure  4b),  which  maximizes  the  reduction  in 
average  MK  height  at  T2.  From  equations  (10)  and  (12),  the  availability  of  mobile 
dislocations  at  the  initial  stage  of  plastic  flow  at  T2  tends  to  be  more  limited  by 
decreasing  /ha  from  0.1  to  0,  and  this  restricted  source  effect,  along  with  the 
diminished  temperature  dependence  of  ha^  due  to  zero  /ha  (equation  (15)),  is 
believed  to  lead  to  the  almost  temperature-independent  0.2%  flow  stresses  in  regime 
A  in  figure  4b  (particularly  for  10%  pre-strain).  Here,  it  was  found  that  equation  (12) 
played  a  crucial  role  in  determining  the  critical  temperature  Tc\  H0  <  h„  (hence 
h0(j,)  —  H0)  for  regime  A  and  hQ  <  H0  (hence  h„{P)  —  h0)  for  regime  B. 

These  effects  can  be  qualitatively  described  as  follows.  In  regime  A,  the  influence 
of  the  built-in  7)  dislocation  substructure  is  so  strong  that  it  remains  in  the 
subsequent  T2  re-straining  and  significantly  diminishes  the  temperature-dependent 
flow  (anomaly)  at  T2.  For  this  case,  the  strength  of  the  stored  substructure  from  7) 
masks  the  flow  anomaly  in  this  regime.  Without  such  a  7)  pre-strain,  the  flow 
anomaly  would  dominate  the  regime.  In  regime  B,  however,  T2  is  high  enough 
that  the  KWL-formation  rate  is  high  and  exhaustion  of  mobile  density  dominates 
flow;  hence,  the  positive  temperature-dependence  of  flow  is  recovered.  For  such 
a  case,  the  stress-controlling  MK  distribution  (mobile  density)  is  dominated  by  the 
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Figure  4.  0.2%-offset  shear  stress  vs.  temperature  plots,  simulated  for  cases  pre-strained  by 
5  or  10%  axial  strain  at  300  K  followed  by  re-straining  at  temperatures  from  300  to  900  K.  (a) 
/ha  =  0. 1  and  (b) /ha  =  0,  where  applicable.  Values  of /hc  and/9c.  were  set  to  0.2  and  0.4  for  5% 
pre-strain  and  0.45  and  0.5  for  10%  pre-strain,  respectively.  The  simulated  original  0.2%- 
offset  shear  stress  versus  temperature  plot  is  also  shown  for  comparison. 


rapid  locking  of  screw-character  dislocations  at  T2  rather  than  by  the  Tx 
substructure. 


Other  comparisons  and  insights 

To  gain  further  insight,  the  simulated  flow  curves  for  selected  temperatures 
in  regimes  A  and  B  were  plotted  together.  Figure  5  shows  simulated  5% 
pre-strained  and  original  flow  curves  at  300  and  900  K,  which  were  used  to  obtain 
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Figure  5.  Simulated  axial  stress  v.s.  axial  strain  curves  re-strained  at  300  and  900  K, 
respectively,  after  being  pre-strained  by  5%  axial  strain  at  300  K.  Values  of /hc,/9c.  and/ha  were 
set  to  0.2,  0.4  and  0.1,  respectively,  for  the  current  simulations,  where  applicable.  Original 
axial  stress  vs.  axial  strain  curves  simulated  at  300  and  900  K  are  also  shown  for  comparison. 


the  r 0.2%-T  plot  in  figure  4a  (/"ha  =  0.1).  The  simulated  5%  pre-strained  flow  curve  at 
300  K  exhibited  a  monotonic  rise  over  an  initial  part  of  the  flow  curve,  with  no 
noticeable  change  in  strain-hardening  behaviour  when  compared  to  the  simulated 
original  flow  curves  at  the  same  temperature.  However,  the  simulated  5%  pre- 
strained  flow  curve  at  900  K  showed  a  micro-yield  transition  behaviour  in  the  axial 
strain  range  between  about  0.3  and  1.4%,  as  indicated  in  figure  5.  This  micro-yield 
transition  seems  to  be  responsible  for  the  slightly  decreased  pre-strained  0.2%  flow 
stresses,  compared  to  those  from  original  flow  curves,  in  regime  B  in  figure  4a. 

As  previously  mentioned,  in  regime  B,  the  built-in  Tx  substructure  effect  through 
Ha  in  equation  (10)  is  no  longer  expected  because  ha  <  Ha  (hence  ho(ji)  —  h0). 
However,  Hco  <  hco  (hence  hco^  —  Hco)  still  holds,  even  in  this  regime.  Note  that  the 
reduction  in  the  critical  MK  height  (/i^))  due  to  T\  pre-strain  leads  to  an  increase  in 
mobile  MK  density  (pmK(/>))  fr°m  equation  (6).  This  can  be  understood  as  an 
increased  population  of  mobilized  MKs  due  to  the  7)  pre-strain  -details  of  which 
will  be  discussed  later.  Now,  initial  plastic  flow  at  T2  tends  to  rely  heavily  on  the 
variation  of  /hc  and  f8c  (equations  (9)  and  (11)).  In  particular,  it  was  discussed  in  the 
previous  section  that  the  micro-yield  transition  range  at  T2  is  controlled  by  f9c 
(figure  2).  Figure  6  shows  the  simulated  5%  pre-strained  flow  curves  at  800  K  for  two 
different  values  of  f8c,  0.4  and  0.6,  along  with  the  simulated  original  flow  curve. 
Initiation  of  the  micro-yield  took  place  at  the  same  axial  strain  for  both  5%  pre- 
strained  flow  curves  having  different  values  of f8c.  However,  a  more  pronounced  and 
extended  micro-yield  transition  was  observed  as  f8c  decreased  from  0.6  to  0.4.  Also, 
upon  completion  of  the  micro-yield  transition  the  simulated  5%  pre-strained  flow 
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Figure  6.  Simulated  axial  stress  v.s.  axial  strain  curves  for  two  different  values  of fgc,  0.4  and 
0.6,  re-strained  at  800  K  after  being  pre-strained  by  5%  axial  strain  at  300  K.  Values  of/hc  and 
/ha  were  set  to  0.2  and  0.1,  respectively,  for  the  current  simulation.  The  original  axial  stress  vs. 
axial  strain  curve  simulated  at  800  K  is  also  shown  for  comparison.  The  insert  is  the  magnified 
micro-yield  transition  regime. 


curves  returned  to  follow  the  original  strain-hardening  behaviour.  Note  that  the 
micro-yield  transition  shown  in  the  simulated  pre-strained  flow  curves  at  900  K 
(figure  5)  and  800  K  (figure  6)  is  in  qualitative  agreement  with  experimental 
observations  for  low  temperature  pre-straining  followed  by  high  temperature 
re-straining  [5,  13]. 

For  the  case  of  Tt  <  73.  the  two-step  deformation  mechanism  interpreted  from 
the  current  simulation  results  (figures  4-6)  was  compared  to  mechanisms  proposed 
by  Ezz  and  FLirsch  [12]  and  Shi  et  al.  [13].  In  regime  A,  Ezz  and  Hirsch  [12]  ascribed 
the  increased  pre-strained  flow  stresses  to  the  mobile  dislocation  sources  that  were 
already  sufficiently  hardened  during  pre-straining  at  Tu  leading  to  ryp(T2)  larger 
than  rs(T2)  in  [12],  where  r/73 )  is  the  athermal  bypass  stress  proposed  by  Hirsch 
[26,  27],  which  was  also  expressed  as  a  function  of  exp(-H//3kT)  [15],  where  77/  is  the 
activation  energy  for  cross  slip.  This  is  probably  in  accordance  with  the  built-in  Tx 
limited-source  effect,  which  was  implemented  in  re-straining  at  73  through  the  initial 
average  MK  height  at  73  (h0(Pj)  in  equation  (12).  This  occurred  because  the 
MK-height  configuration  (hence  the  dislocation  substructure)  introduced  by  pre¬ 
straining  at  7)  (represented  by  II „)  is  stronger  than  that  at  73  without  pre-strain 
(represented  by  h0),  for  which  in  turn  H0  <  ha  (hence  li0(P)  —  Ha).  In  regime  B, 
however,  they  suggested  that  ryp(T2 )  is  almost  comparable  to  rs.(73),  which  implies 
that  the  reversible  nature  of  ts(T2)  dominates  the  pre-strain  effect.  In  the  present 
simulation,  li„{P)  returns  to  the  original  h0  in  regime  B  because  77(,  is  no  longer  strong 
enough  to  defeat  ha  in  this  regime.  Hence,  h0(p)  equals  h„.  This  means  that  it  takes 
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back  the  thermally  reversible  nature  of  ha  (equation  (11)  in  [17]),  which  led  to  the 
recovered  positive  temperature  dependence  of  the  flow  stress  in  regime  B.  They  also 
suggested  that  the  slightly  decreased  r yp(Tf),  compared  to  r y(T2),  in  regime  B  was 
attributed  to  the  propagation  of  long  edge-character  dislocations,  which  was 
facilitated  by  forest  obstacles  inherited  from  pre-straining  at  T\  [12]. 

Furthermore,  Shi  et  al.  [13]  suggested  that  the  micro-yield  transition  in  the  early 
stage  of  plastic  flow  at  T2  was  the  result  of  an  additional  population  of  large  MKs 
built  during  pre-straining  at  Tl .  In  the  present  simulation,  the  micro-yield  transition 
and,  hence,  the  corresponding  flow-stress  drop  in  regime  B  arose  because  the  critical 
MK  height,  determined  from  the  T ,  pre-strain  (represented  by  //,„),  was  still  smaller 
than  that  at  T2  without  pre-strain  (represented  by  hco)  in  this  regime;  hence 
hcotp)  —  Hco.  In  other  words,  some  population  of  larger  MKs  was  carried  from  T \  to 
T2  through  hco(p). 

Although  there  were  slight  disparities  in  the  detailed  understanding  of 
responsible  mechanisms,  interpretations  of  pre-straining  effects  discussed  in  the 
literature  were  in  overall  agreement  with  the  present  simulations.  This  is  particularly 
true  with  respect  to  the  idea  of  a  temporary  ‘readiness’  for  mobilization  and 
propagation  of  edge-character  dislocations,  mainly  MKs  at  the  early  stage  of  re¬ 
straining  at  T2,  being  responsible  for  the  slight  flow-stress  drop  and  the  occurrence  of 
the  micro-yield  transition  in  regime  B.  Also  note  that  there  is  excellent  overall  accord 
between  the  controlling  mechanisms  deduced  from  the  present  modelling  and  those 
previously  proposed  from  experimental  studies  [12,  13].  However,  even  with  such  a 
model  and  clear  evidence  of  MK-controlled  stress,  there  is  no  new  light  to  be  shed 
upon  the  exact  mechanisms  by  which  a  sufficient  population  of  MKs  persists  in  the 
substructure,  even  after  saturation  of  straining  under  sustained  load,  as  performed 
by  Ezz  and  Hirsch  [15].  One  is  only  left  to  conclude  that  some  substructure  relaxation 
always  takes  place  on  unloading  from  Tu  leaving  mobile  segments  in  the 
substructure.  Upon  re-loading  at  T2,  the  stress  to  activate  and  multiply  dislocations 
from  those  segments  over  a  range  of  strain  will  be  controlled  by  either  the  details  of 
the  T |  structure  or  the  KWL-controlled  mobile  density  at  T2,  depending  upon  the 
values  of  Tu  T2,  strains  at  each  temperature  or  their  sequence  in  time. 


6.  Concluding  remarks 

We  have  shown  that  the  two-step  (T\  and  Tf)  deformation  behaviour  of  Ni3Al-based 
single  crystals  can  be  qualitatively  represented  under  the  framework  of  a  new 
constitutive  model  proposed  in  [17].  Pre-strain  at  T \  gave  rise  to  substantial  changes 
of  plastic  flow  at  T2,  such  as  a  significant  increase  in  flow  stress  (At),  weakened 
positive  temperature-dependence  of  flow  stress  and  the  occurrence  of  a  micro-yield 
transition.  All  these  pre-strain  effects  were  found  to  be  qualitatively  but  self- 
consistently  resolved  in  modelling  through  a  reasonable  partitioning  of  reversible 
and  irreversible  contributions  of  the  dislocation  substructures,  i.e.  by  understanding 
the  dislocation  substructures  evolving  at  T\  and  through  controlling  the  strain 
duration  of  these  contributions  at  T2.  In  particular,  a  large  flow-stress  increment  at 
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T2,  despite  only  a  small  magnitude  of  pre-strain  at  7j,  was  achievable  by  controlling 
the  availability  of  mobile  dislocations  and  sources  at  T2. 

The  present  constitutive  modelling  was  able  to  qualitatively  capture  all  such 
behaviours,  since  it  was  able  to  handle  the  density  of  mobile  dislocations,  in 
particular,  a  mobile  MK  density  (pmK),  as  a  direct  constituent  of  plastic  flow.  This 
approach  is  significantly  distinct  from  conventional  crystallographic  constitutive 
models  in  which  an  abundant  mobile  density  is  inferred  and  stress  is  linked  to  the 
Taylor-type  homogenization  (x strength  ~  a  up'12)  rule.  In  those  models,  the  total 
dislocation  density  is  stress-controlling  and  evolves  according  to  phenomenological 
rules  for  storage  and  recovery.  In  the  new  model,  the  mobile  dislocation  density  is  a 
degree  of  freedom  itself,  which  is  directly  linked  to  evolution  statistics  of  the  MK 
height  distribution  and  may  limit  stress,  depending  on  temperature  and  strain 
history. 
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